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A computer program was written for use in sizing gas turbine counter- 
flow regenerators. The program simulates the semi-open gas turbine re- 
generative cycle under constant load conditions. Inputs to the program 
include compression ratio, turbine inlet temperature, component efficien- 
cies, horsepower output and regenerator core geometry. External re- 
generator core dimensions are varied and various outputs are computed 
such as thermal efficiency, mass flow rate, total relative pressure 
drop, regenerator effectiveness and temperatures and pressures at all 



key points in the cycle. 
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1 . 



Introduction. 



The design of a heat exchanger involves a consideration of both 
the heat transfer rates between the fluids and the mechanical pumping 
power expended to overcome fluid friction and move the fluids through 
the heat exchanger. Kays and London state that for a heat exchanger 
operating with high density fluids, the friction power expenditure is 
generally small relative to the heat transfer rate, with the result 
that friction power expenditure is seldom of controlling influence (1). 
However, for low-density fluids, such as the gases flowing through a 
gas turbine regenerator, it is very easy to expend as much mechanical 
energy in overcoming friction as is transferred as heat. The regenera- 
tor designer must therefore establish the proper balance between heat 
transfer rate and pressure drop. 

Two advantages of the gas turbine are compactness and low weight 
whereas heat exchangers are inherently rather bulky and heavy. Weight 
and size considerations must therefore be taken into account along with 
heat transfer rate and pressure drop balance if regenerators are to 
serve a useful purpose. Since the prime reason for regenerators is 
the improvement of overall thermal efficiency, it appeared that a fruit- 
ful approach to regenerator sizing and core selection would be considera 
tion of the gas turbine cycle as a whole rather than the regenerator 
separately. Adjustments in internal regenerator core geometry and ex- 
ternal dimensions are then reflected directly in terms of changes in 
thermal efficiency rather than indirectly in terms of changes in pres- 
sure drop and heat transfer rate. A computer program was written for 
this purpose which essentially simulates the semi-open gas turbine re- 
generative cycle under constant load conditions. Inputs to the program 
are compression ratio, turbine inlet temperature, component efficiencies 
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horsepower output, regenerator core geometry, etc. External regenerator 
core dimensions are then varied and various outputs are computed such as 
thermal efficiency, mass flow rate, total relative pressure drop, re- 
generator effectiveness and temperatures and pressures at all key points 
in the cycle. 

Since the program is based on fundamental relationships from thermo 
dynamics, heat transfer and fluid mechanics, sections 2, 4 and 5 are a 
review of the basic theory utilized. 
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TEMPERATURE, ENTHALPY 



2 . 



Thermodynamics of Semi-Open Regenerative Cycle 




FIGURE 1. SEMI -OPEN REGENERATIVE CYCLE 
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2.1 Description of the Actual and Ideal Cycles. 

Both the ideal and actual regenerative cycles are shown in Figure 
1. The state points indicating the ideal cycle are subscripted "s"; 
the nonsubscripted state points indicate the actual cycle. 

First describing the ideal cycle, air continuously enters the com- 
pressor in condition "Is". It is compressed isentropically to state "2s" 
then enters the cold side of the regenerator. Theoretically, if the re- 
generator is large enough and the flow is slow enough, the air from the 
compressor can be heated reversibly to temperature "5s" at state "3s". 

It next enters the combustor, some of it going around the outside of 
the combustor chamber proper and some furnishing oxygen for burning fuel 
which is continuously injected. Because of the temperature rise, the 
gases expand and enter the turbine in state "4s". After isentropic ex- 
pansion through the turbine, the hot exhaust gases enter the hot side of 
the regenerator. Again, if the regenerator is large enough and the flow 
slow enough the gases are cooled reversibly to temperature "2s" at state 
"6s" and exhausted. 

The ideal cycle can be described as follows; 





(2.1) 




(2.2) 


ItJ J ~[ji "Rts 


(2.3) 



In the actual cycle, air continuously enters the compressor as in 
the ideal cycle; condition "1" is identical to condition "Is". 
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(2.5) 






It is then compressed in an irreversible adiabatic process to state 
arriving at the same pressure as in the ideal case but with increased 
entropy due to fluid friction. 



There will of course be a pressure drop as the air flows through the 
regenerator and the air is heated only to the temperature at state *’3” 
not because less heat is actually transferred than could be trans- 

ferred reversibly. The air and now gases undergo another pressure drop 
in flowing through the combustor j however they arrive at the turbine 
inlet at the same temperature as in the ideal case. 



To continue the analysis it is necessatry to look ahead to state point 
The pressure at this point can be considered constant because it 
is assumed that the turbine is exhausting into a large volume, most 
likely the atmosphere. 



There will of course be a pressure drop through the hot side of the re- 
generator which indicates that must be higher than if P^ is to 
remain the same as in the ideal case. The gases will also experience 
another increase in entropy as they pass through the turbine. Finally, 
since less heat is actually transferred from the hot side of the re- 
generator than could be transferred reversibly, T^ is higher than 




( 2 . 6 ) 




(2.7) 




( 2 . 8 ) 
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2.2 Compressor. 



If the compressor is considered a steady flow machine with kinetic 
energy change equal to zero, the works of the fluid for adiabatic flow 
are: 

ideal: V/cs = (2.9) 

actual: Wc. = ~ jm. (2.10) 



The foregoing equations may be obtained from the steady flow energy 
equation (2). The compressor efficiency is then: 






(L= 



1D£AL 

ACTt/AL FLUID WOfi-K 



k*s ~ hlS 

h-iS — A. 2. 



( 2 . 11 ) 



2.3 



Combustion Process. 




Figure 2 Energy Diagram for Combustor 

The energy diagram shown above for a steady flow combustor, is on 

the basis of one lb air entering and lb fuel/lb air. The 

m /mm 

entering molecular Internal energy and flow work of the air is h^ and 

of the combustible is . The stored chemical energy in the fuel 

q^is measured from the zero datum of enthalpy. Departing energy is the 

enthalpy of (1+ products and that part of the chemical 

energy (1- 7 ^coa»B> that is not converted into molecular 

form. Combustor efficiency as used here is defined as: 

9^ _ actual energy taken on by the gases 

Iw/hD energy released during complete combustion 
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The resulting equation is 



K 3 ^ P (2.12) 

Then 

(2.13) 

(A j- ^7(cof\S>^-i — 

The lower heating value of the fuel is used in gas turbine analysis 

as it is impractical to expand the gas to a point where water vapor from 
the products of combustion will condense and contribute any energy to 
the cycle (3) . 

The enthalpy of liquid octane with respect to a datum of octane 
vapor having h=0 at absolute zero is given by 

A^- (2.14) 

The enthalpy of other hydrocarbon combustibles can safely be taken as 
the same (the amount of energy is small) (1). 





2.4 Turbine. 

If the turbine is considered a steady flow machine with kinetic 
energy change equal to zero, the works of the fluid for adiabatic flow 
are 



ideal: ^kn) <2-15) 

actual: U4 ^(l -V *^^ 4 ) (A 4 _ ( 2 - 1 &) 

As for the compressor, the foregoing equations can be obtained from the 
steady flow energy equation. The turbine efficiency is then 

-y . ACTUAL FLUiO Wotif< _ k^s-k.r /o 

' 1 ^- ideal 
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2,5 Overall Thermal Efficiency. 

The overall thermal efficiency is defined as net work out divided 
by energy chargeable against the cycle. 




(2.18) 



The energy chargeable against the cycle is the energy released by the 
complete combustion of the fuel. 

( 2 . 19 ) 

The equation for thermal efficiency then follows from equations (2.10), 
(2.16), (2.18), and (2.19) 

■yj , _ ( 1 ( 2 . 20 ) 

2.6 Horsepower Output. 

The equation for horsepower output follows directly from equations 
(2.18) and (2.19) 

H P= U4 (2.21) 

where m has units of lb air per hour and 2544 BTU = one horsepower 
a m 

hour. 



2.7 Specific Fuel Consumption. 

Specific fuel consumption is defined as lb fuel burned per unit 

m 

work output. If the horsepower-hour is used as the unit of work then 
specific fuel consumption has the units of Ib^ fuel per horsepower-hour. 
It follows that 



5FC = 



2T44 



( 2 . 22 ) 
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3. Regenerator Surface and Core Geometry. 

Certain geometrical relations are necessary in the application of 
the basic heat transfer and flow-friction data to a design problem. One 
particular set and form of these relations that has proven convenient 
for use in the computer program is given here. 




Figure 3 Simple Plate -Fin Counterflow Regenerator Core 
The following part of basic data are given for all surfaces in Kays 
and London (4) . 

b = plate spacing 

Vh = flow passage hydraulic radius 

^ = ratio of total transfer area of one side of 

exchanger to volume between plates of that side 

a = plate thickness 

In addition, the ratio of fin area to total transfer area is normally 
given. 

The following geometrical factors are required for each of the two 
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sides, hot and cold, of the complete heat exchanger core. 



A total transfer area of one side of exchanger 



A =» free flow area of one side 



c 



Afy~ = frontal area of one side 



L = flow length of one side 



V = total exchanger core volume 



= ratio of total transfer area of one side of exchanger 
to total exchanger volume 

0” = ratio of free-flow to frontal area of one side of 
exchanger 

V = solid core volume 
s 

The equations below give the relations between surface and core 
factors for one side of the exchanger. Subscript "1" refers to any one 
side and "2" refers to the other side; factors without a subscript are 
common to both sides. 




(3.1) 




L Aiv- 



A i-K 





(3.3) 




(3.4) 




(3.5) 




22 



A,= 



(3.7) 



(3.8) 
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4. Counterflow Exchanger Heat Transfer Theory 
4.1 Conduction Heat Transfer. 

When a temperature gradient exists in a body there is an energy 
transfer from the high temperature region to the low temperature region. 
The energy is transferred by conduction and the heat transfer rate per 
unit area is proportional to the normal temperature gradient. 




(4.1) 



When the proportionality constant Is Inserted, 






(4.2) 



where "k" is called the thermal conductivity of the material. 




Figure 4 Elemental Volume for One-Dimensional Heat 
Conduction Analysis 



If the steady state one-dimensional case is considered 




(4.3) 



which integrates to give 



/aCT.-TJ , 

L 



(4.4) 



providing thermal conductivity is constant. 



4.2 Convection Hi ‘ ' 




Figure 5 Sketch Illustrating Convection from a Plate 



24 



The convective heat transfer rate is related to the overall temperature 
difference between the wall and fluid and the surface area as follows: 

A A (Tw ■" To) (4,5) 

where ”h*' is called the convection heat transfer coefficient. It may be 
anticipated that convection heat transfer will have a dependence on the 
thermophysical properties of the fluid (thermal conductivity, specific 
heat, density, viscosity), the speed at which the fluid is moving and 
some characteristic dimension. If a dimensional analysis is performed 
using h; u, Cp, ^ , k^, ^ and L as variables, the variables may be 
arranged into three independent dimensionless parameters. The para- 
meters depend on how the analysis was performed. One very useful set of 
parameters is 



k r P L C 4^ 



( 4 . 6 ) 



A function of the form 

= (9 (^ 7 ) 

must therefore exist. Such a function may be determined experimentally 
in the form 

< 4 - 8 ) 

for various heat transfer surfaces and geometries. The generalized 
grouping Pv* ^ is defined as the Colburn j-factor and the char- 

acteristic dimension "L" in "Re" is designated as the hydraulic dia- 
meter. 
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p ^ (4.9) 

A large number of curves expressing the relationship given by equation 
(4.8) are published in Kays and London (1). One of the more common ex- 
perimental methods for obtaining these characteristics is described in 
a paper by Pucci, Howard and Piersall (4). 

After determining the Colburn j factor, St is easily found pro- 
viding the fluid properties necessary to determine Pr are obtainable. 

W-IO) 

The convection heat transfer coefficient is then evaluated by using 



K = Gfs 



4.3 Fin Effectiveness. 



t; 




(4.11) 



Figure 6 Sketch Illustrating One -Dimensional Conduction 
Through a Rectangular Fin with Convection at 
the Surface 

Consider the one dimensional fin exposed to a surrounding fluid at 

a temperature T oo • The temperature of the base of the fin is T . 

o 

Making an energy balance on an element of the fin of thickness dx as 
shown in figure: 

(energy in left face)=(energy out right face)+(energy lost by convection) 
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. ^ A 

energy in left face = ^ 

^ 'j' \ 

energy out right face = ->^A j 






(4l+e<il^ 

\t(x «<T 



Let P equal the perimeter of the fin then; 

energy lost by convection = ~ -A. P<^X ( l~"'Too'j 

Conbining quantities, the energy balance yields 



zXv>- JA I' ^ 



(4.12) 



(4.13) 



a a 

Let B - T - Tbc . Then (4.12) becomes 

Ax’- iA ^ ° 

One boundary condition is 

e- 

For the case of a fin with an insulated end, dT/dx = 0 and d Q /dx = 
0 at X = Z. 

Let 

xon •= 



AP 



lA 



The general solution for (4,13) may be written 



(4.14) 



_.»nX 

B ~ L\ e. •¥ Ct.^ 



(4.15) 



Solving for constants and the solution is obtained as 



A. 

Qa 






A, 

' e 



+ 



<2. 



/j^X 



\ + e. 






(4.16) 



The total heat lost by the fin must equal the energy conducted through 
the base of the fin 






X -o 
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Then 









L h 



+ iL 



I 'h 



(4.17) 



R/Aa^ 

Fin efficiency is defined as 

- actual heat transferred 

heat which would be transferred if entire fin area were at base 

temperature 



which becomes 



X?'h9o 



(4.18) 



To apply equation (4.18) in the computer program it is necessary to 
put it in terms of the heat exchanger core geometry. Referring to Figure 



6 , 



so that 



p=:XL-\-lS- ^ L 
A - ^ L 



It is seen from Figure 3 that the fins of the heat exchanger are 
symmetrical with respect to a plane midway between the plates separat- 
ing the fluids. Therefore, the temperature gradient in the fin at this 
plane is zero and the above insulated end solution applies. 

Noting that Z ^ b/2, the fin efficiency becomes 



%- 






(4.19) 



4.4 Surface Effectiveness. 

The total effectiveness for a finneci surface may be defined as the 
ratio of the total heat transfer of the combined area of the surface and 
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fins to the heat which would be transferred if this total area were 



maintained at the base temperature Let "A** equal the total heat 

transfer surface including the fin area then 



^0 = 1 - 



(4.20) 



4.5 Overall Coefficient of Heat Transfer. 




Figure 7 Sketch Illustrating-One Dimensional Conduction 
and Convection through a Plane Wall 

Consider the plane wall shown above exposed to a hot fluid on one 
side and a cooler fluid on the other side. The steady state heat trans- 
fer is expressed by 






(4.21) 



The heat transfer process may be represented by the resistance network 
shown and the overall heat transfer is calculated as the ratio of the 
overall temperature difference to the sum of the thermal resistances. 



Th -Tc 

+|;|- 



(4.22) 



The overall heat transfer coefficient "U", is defined by the re- 



lation 



/A A To 



Vis /CALL 



(4.23) 
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where A is some suitable area for heat flow. Then 

'u. + 4 /A ' 

Now consider a plate of thickness "a" with extended surfaces or fin 
areas on either side. The overall heat transfer coefficient becomes 

U= — ^ ^ ^ — 






a 



(4.24) 






(Ac.//Au) 

where U is based on a unit of hot side total area including fin or ex- 
tended area. denotes the average wall area. In the computer program 

the very small wall resistance is neglected in calculations so that 



U 



-f- 



1 



4.6 



A w JfcyC Ac 

Exchanger Heat Transfer Effectiveness 




fiC/T 

TCj /yV 



Figure 8 Temperature Profiles in Counterflow 
Heat Exchanger 

The actual heat transfer in a heat exchanger may be computed by 
calculating either the energy lost by the hot fluid or the energy gained 
by the cold fluid. Referring to Figure 8, 

(4.25) 

Heat exchanger effectiveness is defined as 



e = 



actual heat transfer 



maximum possible heat transfer 
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To determine the maximum possible heat transfer for the exchanger it is 
necessary to recognize that this maximum value could be attained if one 
of the fluids were to undergo a temperature change equal to the maximum 
temperature difference present in the exchanger, ( Xcy /// )• 

The fluid which might undergo this maximum temperature difference is the 
one having the minimum value of since the energy balance re- 

quires that the energy received by one fluid is equal to that given up 

» 

by the other fluid. If the fluid with the larger value of goes 

through the maximum temperature difference, the other fluid would under- 
go a temperature difference greater than maximum. The maximum heat 
transfer is therefore expressed as 













(4.26) 



The ’’minimum’* fluid may be either the hot or cold fluid depending on the 
mass flow rates and specific heats. 



h 

(Tv4^/V - Tt/ /V ) 



(4.27) 



An expression for effectiveness in counterflow will now be derived. Re- 
ferring to Figure 8, the heat transferred through an element of area dA 
may be written 






(4.28) 



The heat transfer could also be expressed as 

- (J ( Th -"Tc)u^ A (4.29) 

Equations (4.28) and (4.29) combine to give 
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^(Th-Tc) ^ 
Th - Tc 



_ 



Ui<A 



Integration yields 



/J (T\AjJ^ ^ _ 

It \ ~ 

V ‘H)6UT- Tc^ /V j 



or 



(TwjOuJ \Cy/A/ ) 

(Th^/A/'“ Tcy<5 0l) 

C\t\^ dc C \A ' C 



^ c 



] C 






u^ 


C c. 


Cn 


. Cu 


1 UA 


c7 


J 


Cc " 


|c,(3Ur—Tc/ /A/ 



Then 



and 



Th/,aJ-Tc^,^ Cw Tw,(?6>t 



\c.jOuT - T<v^'<^ + ^ ^ (^Tw^/-</~Te/W) 



From (4.32) 



or 



Tuy/v - ^ ( Th/ /A/ - Tcy/A/ ) — Tcy //j _ - |j~ Xt\ 



- £ 



- 6 



Cc 



' cTJ 



(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 

(4.35) 

(4.36) 
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Solving for 






Cu~\ ^ 









If Ca^Ci\ or Cc 






In either case 



I - c« 

the same procedure yields 

UA 

-L'-^trJ 

- e 



£ = 



-b- &] !. 

TcTTTT^l ^ 

r u ^ L J Cc 



fi_ I 

- !_' C/nA'xJ C/*i'a/ 



£=• 



I _ „ ( I - 1 UA. 

' d>«AX 2 ^ Ca^axJcIa.-/ 



1 - 2 



r -Ciii^lvTU 

C A 1 AX J 



C rAihJ p 

(C./v)AX <L 



C 



•j/VT 



U 



(4.37) 



(4.38) 



(4.39) 



(4.40) 



where 

is defined as the number of heat transfer units. In general it is 
possible to express 





6 




J Flow arrangement) 
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5. 



Exchanger Flow-friction Theory. 




Figure 9 Heat Exchanger Core Model for Pressure -Drop 
Analysis 



5.1 Core Friction. 

Friction factor is defined on the basis of an equivalent shear 
force in the flow direction per unit of wetted area. For most of the 
surfaces considered it is a combination of viscous shear and pressure 
force but no attempt is made to separate these two effects. The pres- 
sure force (form drag) is most significant in interrupted or wavy fin 
surfaces where fluid separation from these surfaces occurs. Thus 

( 5 . 1 ) 

or 



f = 1 . 3 ^° ^ 



(5.2) 



It follows directly then that 



AP-i 



k P 









A 



(5.3) 



Ac Ac 1 ^ A 

where equals the average density in the core and "G" is based 

on free flow area- Curves showing friction factor as a function of 
Reynolds number are published in Kays and London (1). They are plotted 
for the same surfaces as the Colburn j- factor mentioned in the previous 
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section. 



5.2 Flow Acceleration. 

A term accounting for flow acceleration must also be included be- 
cause of the change in density as the fluid flows through the core. A 
force balance using change in momentum yields 

(5.4) 



or 



Ap- i\ 

^ [p^ Ij 



(5.5) 



The core friction term and the flow acceleration term may now be com- 
bined to give the pressure drop between points "a" and "b". 



f A ^ ±_( Po. 






f 



(5.6) 



5.3 Entrance and Exit Effects. 

A typical gas turbine regenerator installation usually involves a 
flow contraction at the heat exchanger core entrance and a flow expan- 
sion at the core exit. They introduce additional flow stream pressure 
drops. This problem has been investigated both analytically and experi- 
mentally by Kays (5) . 

The entrance pressure drop is made up of two parts. The first, 

(P^ - Pg)j^» is the pressure drop which would occur due to a flow area 
change alone. The second, (P - P )., is the pressure loss due to the 

X 8L 

irreversible free expansion that always follows the abrupt contraction. 

First considering the pressure drop due to area change alone, 

liJr (5.7) 
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The pressure change from "1" to "a" is very small relative to the total 



pressure so that p a . Therefore 



A ~~ /^e. lA-o. 



(5.8) 



Substitution of equation (5.8) into (5.7) yields 

The analysis, as presented in Kays, of the pressure drop due to 

free expansion is somewhat more complex and only a brief description of 

the equations utilized in the computer program will be given here. The 

contraction loss coefficient, K is defined as follows: 

c 



(r-pA 



-K ^ - k, 



6-’ 






(5.10) 



K can be evaluated from 
c 












(5.11) 



is a function of Reynolds number and flow geometry. For laminar 
flow 



= 1.333 (circular tubes) 
K, = 1.390 (square tubes) 



= 1.430 (triangular tubes) 



(5.12) 



For turbulent flow in circular tubes 






(5.13) 



where 



0.04'i Re 



-0.2 



(5.14) 



and 



kU(SWAlie') -- I _i iq 

KaL(Cllie.ULAl^)~\ ^ 



(5.15) 
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is a function of the contraction ratio CT and is presented graphical- 
ly for circular tubes and a one-dimensional contraction. Kays recommends 
that the one-dimensional contraction be used in evaluating triangular 
and square tubes. The curves of vs. CT were fit with straight line 
approximations for use in the computer program and are shown in Ap- 
pendix B. 

The exit pressure rise is similarly made up of two parts. The 
first, (P- - Po)-i> is the pressure rise which would occur due to area 
change alone. The second, (P^ - 1^2^ 2^ pressure loss associated 

with the irreversible free expansion and momentum changes following an 
abrupt expansion. The first term follows from the same analysis as the 
corresponding term in the entrance pressure drop. 

(P^- PA = ^ A"-') 

The free expansion term is given by 

T_ 

(p^ - R.) x _ ^ (5.17) 

where 

\-ZKd(T^ r"' (5.18) 

and K, is evaluated using (5.12), (5.13) and (5.15). 
o 

and are presented graphically in Kays and London for the 
circular, triangular and square tubes; however, the preceding equations 
were used for evaluation in the computer program. 

If the core friction, flow acceleration, and the entrance and exit 
effects are now combined 
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h- ^ 

Ac /^ 




Vi 

5-t 



(l-<r'-+/f,') + ^(% 




+ f A 

Ac^i ^ 



(o-“l + 




(5.19) 
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6. Computer Program Development. 

6.1 General Approach. 

As discussed in the introduction, the program simulates a gas tur- 
bine arrangement running under constant load conditions ♦ The nature of 
the problem necessitates that an iterative approach be taken in solving 
for most of the desired outputs. The program is flexible in that six 
different sets of output data are obtainable depending on what type of 
information the user desires. 

Program inputs can be divided into four categories as follows: 

a. known information describing the cycle as a whole 

b. information describing regenerator core characteristics 

c. assumed quantities necessary to start the solution 

d. special inputs depending on the type of output desired 
The first three categories are referred to hereafter as general inputs. 

Program outputs can be divided into two categories as follows: 

a. special outputs depending on which set of output data is 
desired 

b. general outputs common to all sets of output data 

The program is executed largely through the use of subroutines so 
that characteristics of various gas turbine arrangements are easily 
Interchangeable. All general outputs are solved for in one rather 
long subroutine called XMAIN. Subroutine XMAIN utilizes all the other 
subroutines as the solution progresses. Subroutine XMAIN is described 
in detail later in this section; however, the general approach for the 
solution of the various unknowns in the gas turbine cycle is as fol- 
lows : 

a. Referring to Figure 1, the temperature and pressure at 
point are solved for knowing the temperature and pressure at point 
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**1", the compression ratio and the compressor efficiency. 

b. The pressures at points **5** and **3** are computed using 
assumed values for pressure drops through the regenerator. 

c. An estimate of the temperature at point **5** is computed 



the estimated temperature at point **5** and an assumed value for re- 
generator effectiveness. 

e. Better values are computed for the regenerator pressure 
drops using theory from section 5 and an assumed mass flow rate. 

f. Better values are computed for the temperatures at points 



g. A new value is computed for regenerator effectiveness 
using the theory from section 4 and fed back into step **d**. 

h. The process is repeated until all temperatures and pres- 
sures converge. 

i. A new value of mass flow rate is now computed by comparing 
the calculated value of power output to the desired power output. 

j. The entire process is then repeated until mass flow rate 

converges . 

k. All remaining general outputs are now computed. 

6.2 General Inputs. 

a. The following inputs describe the cycle as a whole: 



using 




and the compression ratio. 



d. The temperatures at points **3** and **6** are computed using 



’’3** and 



T^ compressor inlet temperature 

compressor inlet pressure 
compressor pressure ratio 




C compressor efficiency 
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^4 


turbine inlet temperature 


^6 


hot side regenerator exhaust pressure 




turbine efficiency 


4 As P 

P 3 

HP 


combustor efficiency 

relative pressure drop through combustor 
net horsepower output 


^f 


fuel temperature at combustor inlet 


‘lx 


lower heating value of fuel 


R 


gas constant for air 



b. The following inputs describe the regenerator core character 
istics: 



P 


density of regenerator material 


4 


cold side heat transfer area/volume between 
cold side plates 




hot side heat transfer area/volume between 
hot side plates 


k 


cold side plate spacing 


u 


hot side plate spacing 




cold side hydraulic radius 




hot side hydraulic radius 


Sc, 


cold side fin thickness 




hot side fin thickness 


Af,^ 


cold side fin or extended area/total cold side 
transfer area 




hot side fin or extended area/total hot side 
transfer area 


a 


plate thickness between cold and hot sides 


k 


fin thermal conductivity 



c. The following are assumed input quantities: 

m mass flow rate of air 

a 
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relative pressure drop through cold side of re- 
generator 

relative pressure drop through hot side of re- 
generator 

regenerator effectiveness 

6.3 Subroutines Other Than Subroutine XMAIN 

a. Subroutine KEKC (SIGMA, RE, XKE, XKC) 

Input arguments are CT and Re ; computes Ke ' and Kc 
using equations from section 5. 

b. Function PRAND (T) 

Input argument is T ; computes Prandtl number of air from 
equation obtained by fitting curve to data given in Hilsenrath (6). 

See Appendix B for data. 

c. Function VIS (T) 

Input argument is T ; computes of air from equation 

obtained by fitting curve to data given in Hilsenrath (6). See Ap- 
pendix B for data. 

d. * Function SPEC (A) 

Input argument is T ; computes c of air. The equations 

P 

used for computation in this and the next four subprograms are given in 
Howard (2). 

e. Function TEMR (B) 

Input argument is Pr ; computes T of air. Pr as used here 
refers to relative pressure, not to be confused with Prandtl number. 

f. Function ENTH (B) 

Input argument is T ; computes h of air. 

g. Function TEMH (B) 

Input argument is h; computes T of air. 



3 A2.K 

Pv. 

6ArP 

Pr 

6 
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h. Function PREP (B) 

Input argument is T; computes Pr of air. 

i. Function FUENTH (T) 

Input argument is T^; computes enthalpy of fuel from equation 
given in section 2. 

j . Function FRICC (RE) 

Input argument is regenerator cold side Re; computes cold 
side f from equation obtained by fitting curve to experimental data. 

See Appendix B for details on f and j vs. Re data used in this and the 
following three subprograms. 

k. Function FRICH (RE) 

Input argument is regenerator hot side Re; computes hot side 
f from equation obtained by fitting curve to experimental data. 

l. Function STPRC (RE) 

Input argument is regenerator cold side Re; computes cold side 
j from equation obtained by fitting curve to experimental data. 

m. Function STPEUI (RE) 

Input argument is regenerator hot side Re; computes hot side 
j fran equation obtained by fitting curve to experimental data. 

n. Function GENFX (TU, RATIO) 

Input arguments are C . /C and NTU; computes ^ from 

mrn max 

equation developed in section 4. 

o. Subroutine XWRITE 

This subroutine contains all general output print statements. 

p. Subroutine XOPEN 
Computes open cycle outputs 



43 



6.4 General Outputs. 

r 

m mass flow rate of air 

a 

^ regenerator effectiveness 

fuel to air ratio 

SFC specific fuel consumption 

V regenerator core volume 

Wt. regenerator core weight 

compressor discharge pressure 
T2 compressor discharge temperature 

cold side regenerator outlet pressure 
cold side regenerator outlet temperature 
turbine inlet pressure 
P^ turbine exhaust pressure 

turbine exhaust temperature 
Tg regenerator hot side outlet temperature 

Re^ regenerator cold side Reynolds number 

Re^ regenerator hot side Reynolds number 

NTU number of exchanger heat transfer units 

Cmin/Cmax capacity-rate ratio 



aP \ — 3A2. P , 4A3 P , ^ArP 

P /total R 



6.5 Detailed Developmentof Subroutine XMAIN 

Assumptions made in the development of subroutine XMAIN are: 
a. Combustion products have a negligible effect on the 



fluid properties of gases leaving the combustor. 

b. % . ^ ^ " remain constant dur- 

'3 

ing the iteration process. 
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and 



The first 

and 

and . 


step in the subroutine is to calculate T 2 


from Tj^, Pj^ 




(6.1) 






(6.2) 






(6.3) 




Ts"" M fni!2s) 


(6.4) 






(6.5) 




hi Ai 


(6.6) 




T.= ^(A.) 


(6.7) 


Next, h, and P 
4 


^are computed from T, . 
rel^ 4 








(6.8) 




^^^4 ~ ^ ( T 4 ) 


(6.9) 


and P^ are 

Ar P 


now computed using the assumed values for 


3Az P j 

Pi 


Pr 


R - R( 


(6.10) 






(6.11) 






(6.12) 



It next becomes necessary to get an estimate of T^. This is done by using 
a cycle with no regenerator and no pressure drop through the combustor. 
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(6.13) 



Si 

f? 



s. 

Pr. 




(A.4- 






(6.14) 

(6.15) 

(6.16) 

(6.17) 

(6.18) 



and Tg are now computed using the value just computed for and 
and assumed value of ^ 

"^3 ^ ^(Tr^Tl) ”^Ti. (6.19) 

T< = Tr - i (Tj— T-3 (6.20) 

Temperatures and pressures at points *’2”, **3”, **5" and ’’6" are now 
known or have been estimated. Better values are computed for the pres- 
sure drops through the regenerator using the theory from section 5, an 

assumed value of m and fluid properties computed from the temperatures 

a 

and pressures at the above four points. For the regenerator cold side: 

( 6 . 21 ) 

( 6 - 22 ) 

— ( 6 . 23 ) 

.y.y.±.y±- ( 6 . 24 ) 

■2- 
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V ~ La As(r 


(6.25) 


o(^- 

*4" // A. +“ 1 . 0 . 


(6.26) 


Ac = «^cXV 


(6.27) 


(Tc. -= 0(cX 


(6.28) 


> 

Ji 

X 

>> 


(6.29) 


Gc = ^VAc,. 


(6.30) 


/^C ^('^yc) 


(6.31) 


0 4 (S-c 


(6.32) 


ft = 


(6.33) 



R-f?- 



6-c. 






(\-(Jc -V- ^ ( 




% 





(6.34) 



The pressure drop through the regenerator hot side is 
same manner except that the fuel to air ratio must be 



culating G^. 



/- _ (^l "V 

At 



computed in the 
included in cal- 

(6.35) 



where 




f/i4~ A 3 ) 
(Kj-+ 



(6.36) 
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and are now recalculated using the results just obtained. 



R-(Pa-R') 


(6.37) 


R-R+(Pr-B") 


(6.38) 


A new value of is computed as follows: 




R- 


(6.39) 




(6.40) 


Trs- 


(6.41) 


^ 4 ("n"s) 


(6.42) 




(6.43) 


Tr “ ( ^f) 


(6.44) 



Using this value of and the assumed 6 , new values for and Tg 

are computed. \,h. ’ Mk » ^ 3 , Y i/cK > 

are then recomputed using the above values of and 

J o 

A better £ can now be determined using the theory from section 



?- 

II 


(6.45) 




(6.46) 




(6.47) 


R-_, c - {J'^j c-') 


(6.48) 




(6.49) 
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(6.50) 







tU{\lU/A%ii Hr^ 

fTWTi: 

'yi _ 7^-v«^ ( kc/.£ Sc.' -^c 

^ ^ 

! 1 

Mtk, 



(6.51) 

(6.52) 

(6.53) 

(6.54) 

(6.55) 

(6.56) 

(6.57) 

(6.58) 

(6.59) 

(6.60) 
(6.61) 
(6.62) 
(6.63) 
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The above value of £ is now fed back into equation (6.19) and the 
entire set of calculations from (6.19) through (6.63) is repeated until 
the desired degree of convergence is obtained. Many different criteria 
for determining convergence could be used; however, numerous trials in- 
dicate that the regenerator cold side pressure drop converges at a slower 
rate than any of the other computed outputs. When the value of (P^ - ^ 3 ) > 
computed from equation (6.34), converges to within one percent of its 
previous value the subroutine leaves the above loop and proceeds to the 
next step. 

Thermal efficiency and the true value of mass flow rate needed to 
give the desired horsepower output at that thermal efficiency are now 
computed. 

Then 






(6.65) 



I 

where rn^is the assumed value. The above calculated value of horsepower 
is then compared with the desired value. If the two agree within .1 
percent of the desired horsepower the subroutine ends; if they do not a 

0 

new value of m is calculated 
a 



z 4 HP 



( 6 . 66 ) 



where HP is the desired horsepower output. This value of m^ is then 
fed back into the subroutine where the process of calculating pressure 
drops, effectiveness, thermal efficiency and horsepower is repeated 
until horsepower converges. Convergence is normally rapid and two or 
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at most three times through the set of calculations is sufficient. 

The subroutine then proceeds to calculate the remainder of the 
general outputs. 



SFC - 

Vs =• L [Af V-- Ac.d'iJ 



^ /^XVr 



aP\ 

p L:. ft 






(6.67) 

( 6 . 68 ) 

(6.69) 

(6.70) 



6.6 Special Inputs and Outputs. 

As mentioned before, the program is flexible in that different sets 
of output data are obtainable depending on what type of information the 
user desires. The five possible sets of output and corresponding inputs 
are described below. 

a. Output type one 

In some cases the frontal area and regenerator core 
length may be fixed at some predetermined values. The user may then 
determine gas turbine performance for that frontal area 

and length using output type one. The frontal area and length along with 
the general inputs are read into subroutine XMAIN and all general out- 
puts are computed. 

Inputs: general, length, frontal area 

Outputs: general 

b. Output type two 

It may be desirable to fix the overall thermal efficiency 
at a predetermined value. There may also be some latitude possible in 
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choosing the frontal area of the regenerator providing the length re- 
mains between specified limits. Output type two is used to determine 
a frontal area for each length under consideration that will give the 
desired overall thermal efficiency. The shortest length under considera- 
tion, the longest length under consideration, the size of the increments 
desired between the shortest and longest lengths, the desired thermal ef- 
ficiency and all general inputs are read into the computer program. Sub- 
routine XMAIN is then called with the shortest length and an initial 
value of frontal area, determined by the computer, as arguments. Frontal 
area is then incremented with length held constant until the computed 
thermal efficiency is in agreement with the desired thermal efficiency. 

The same procedure is repeated for all other lengths under considera- 
tion. 

Inputs: general, shortest length under consideration, 

longest length under consideration, size of increments desired between 
shortest and longest lengths, desired thermal efficiency . 

Outputs: general, required frontal area for each length 

considered . 

c. Output type three. 

Output type three is similar to two in that the overall 
thermal efficiency is set at a predetermined value. In this case, how- 
ever, it is assumed that there is some latitude possible in choosing the 
length providing the frontal area remains between specified limits. A 
length must therefore be determined for each frontal area under considera- 
tion that will give the desired thermal efficiency. The same computing 
procedure is followed as in output type two except that frontal area is 
held constant and length is incremented until the desired thermal effi- 
ciency is obtained. 
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Inputs: general, smallest frontal area under considera- 

tion^ largest frontal area under consideration, size of increments de- 
sired between smallest frontal area and largest frontal area, desired 
thermal efficiency. 

Outputs: general, required length for each frontal area 

considered. 

d. Output type four. 

In some cases it may be desirable to find the minimum 
frontal area that will allow a predetermined thermal efficiency. In 
other words make the frontal area as small as possible without giving 
consideration to length. Output type four is used to determine this 
minimum frontal area and the length that must be used with it to obtain 
the desired thermal efficiency. 

All general inputs along with the desired thermal efficiency are 
read into the program. Subrouting XMAIN is called and an initial length, 
depending on the mass flow rate and internal core geometry, is deter- 
mined by the computer. Frontal area is then incremented until the de- 
sired thermal efficiency is obtained. A new length is then selected and 
again frontal area is Incremented until the frontal area corresponding to 
the new length gives the desired thermal efficiency. This process is 
repeated and frontal areas are compared until a minimum allowable front- 
al area is obtained. 

Inputs: general, desired thermal efficiency 

Outputs: general, minimum frontal area for desired 
thermal efficiency, length to be used with that frontal area 

e. Output type five. 

The user may find it helpful to determine the maximum 
thermal efficiency possible with a fixed allowable frontal area or range 
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of frontal areas. Output type five computes the maximum thermal ef- 
ficiency and the length that must be used with the frontal area under 
consideration to obtain that thermal efficiency. The frontal area under 
consideration is read into subroutine XMAIN and length is incremented 
until thermal efficiency reaches a maximum. 

Inputs: general, smallest frontal area under considera- 

tion, largest frontal area under consideration, size of increments de- 
sired between smallest frontal area and largest frontal area. 

Outputs: general, length corresponding to each frontal 
area considered for maximum thermal efficiency, 
f. Output type six. 

Output type six is utilized in determining thermal ef- 
ficiency, mass flow rate, specific fuel consumption, fuel to air ratio, 
compressor discharge temperature and pressure and turbine exhaust tempera- 
ture of the simple open cycle. This information can be used in compari- 
son with other output types in setting a standard for regenerator per- 
formance . 



Inputs: general except regenerator characteristics 
Outputs: as given above 

g. Tabular summary of special inputs and outputs 



Output 

Type 


One 


Two 


Three 


Four 


Five 


Six 


A, 


I 


» « « 
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I 




fr,rain 

A. 
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I 




fr,max 

L 
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0 
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0 








min 
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0 








max 
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0 




opt 


0 
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I 
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0 
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— 



54 




FIGURE 10. FLOW CHART OF COMPUTING PROCEDURE 
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7. 



Results ■ 



The ORENDA OT-4 gas turbine was utilized as a model in the develop- 
ment and testing of the computer program. This engine meets a require- 
ment that was generated by the U. S. Army and Navy in 1960 and is dis- 
cussed in a paper by Quan (7). The design target at 600 hp on specific 
fuel consumption was .40 Ibm of fuel per hp hr, which is equivalent to 
a thermal efficiency of 34.6%. The pressure ratio ranges from four to 
five. The following were additional design targets: 



compressor efficiency 


86% 


turbine inlet temperature 


2200°R 


turbine efficiency 


89% 


combustor efficiency 


98% 


regenerator effectiveness 


82.5% 


A few additional quantities were required to 


complete the general inputs 


in the computer program. Their values were 


selected as realistically as 


possible as follows: 




compressor inlet temperature 


560®R 


compressor inlet pressure 


2120 Ibf-ft"^ 


fuel inlet temperature 


560“R 


lower heating value of fuel 


19,000 BTU-lbm"^ 


regenerator hot side exhaust 
pressure 


2120 Ibf-ft"^ 



Six cores were selected at random for evaluation and incorporated 
into the above gas turbine arrangement. The characteristics of the 
first five shown on the following pages are from Kays and London (1); 
the last one from Pucci et.al.(8). The friction factor and Colburn 
j-factor data utilized in obtaining functional relationships for each 
core is listed in Appendix B along with the curve fitting results. The 
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cores were made of stainless steel, aluminum or nickel; however, the 

density of stainless steel was used as the input in the program for all 

cores to obtain a common basis for weight comparison. Output types one 

through six were computed for each core. The desired thermal efficiency 

of .34 was set for output types two, three and four. Required frontal 

areas were computed to match lengths of 3.0 ft through 7.0 ft in output 

type two and required lengths were computed to match frontal areas of 
2 2 

3,0 ft through 7.0 ft in output type three. All core incrementing in 
output type three was done in ,1 ft increments so that lengths determined 
for = .34 are only accurate to within .1 ft. The thermal efficiency 

obtainable with these computed lengths may therefore be slightly higher 
or lower than .34, The characteristics of each core along with results 
selected from the various program output types are tabulated. 

The tabular results are followed by six sets of curves which il- 
lustrate the range obtainable in overall dimensions and weight depend- 
ing on the core selected. The set of curves shown in Figure 12 illus- 
trates the lengths that must be used over a range of frontal areas to 
obtain a thermal efficiency of ,34. Figure 13 gives regenerator core 
weights as a function of frontal area to obtain a thermal efficiency of 
.34 and Figure 14 is a plot of regenerator effectiveness vs. total re- 
lative pressure drop for a thermal efficiency of .34. Figure 15 il- 
lustrates the maximum thermal efficiency obtainable as a function of 
frontal area and the optimum lengths for each frontal area, l.e,, the 
lengths necessary to obtain maximum thermal efficiency are shown in 
Figure 16. Core weights as a function of frontal area for maximum 
thermal efficiency are shown in Figure 17. 

A sample of the six computer output types, generated using core 1, 
is shown following the curves. Output type one is shown on page 73, 
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two on pages 74 through 76 , three on pages 77 through 79 , four on 
page 80 , five on pages 81. through 83 and six on page 84 . 

As mentioned above, the six cores evaluated were selected at ran- 
dom. No systematic attempt was made in finding an optimum core from 
the standpoint of external dimension or weight considerations. The 
results serve only to illustrate the usefulness of the computer program 
in selecting the best core for a given gas turbine. A few very general 
observations regarding the geometries investigated are as follows: 

1. The desired thermal efficiency was not obtainable with the 
fine mesh cores in the small frontal area range, regardless of length. 

2. Significantly longer lengths were required to produce 

the desired thermal efficiency when the coarse mesh cores were utilized. 

3. The coarse mesh cores add considerably more weight and 
volume to the gas turbine arrangement than the fine mesh cores for a 
given thermal efficiency. 

4. Figure 12 shows that required core weight decreases or 
remains essentially constant as frontal area is increased with 7ft A- " *34 
for all cores except core 4. It is apparent that this should be the 
case after inspection of Figure 11 which reveals that required volume 
remains essentially constant as frontal area is increased for all cores 
except core 4. 

5. Total relative pressure drop vs. regenerator effective- 
ness appears to be a linear relationship for fixed thermal efficiency 
Independent of internal core geometry. 

6. Fine mesh cores are much more sensitive to changes in 
external dimensions than coarse mesh cores. Shorter length Increments 
than .1 ft are necessary for good accuracy when investigating fine mesh 
cores. This is apparent from Figure 13 where the points determined for 
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cores 3 and 6 are the most seriously in error. Instructions for changing 
the length increments in the computer program are given in Appendix A. 

7. Higher maximum thermal efficiencies are obtainable 
with the coarse mesh cores, however, the corresponding lengths and weights 
are much greater than with the fine mesh cores. 

8. Cores 3 and 6, having approximately the same plate spacing, 
produced essentially the same results even though core 6 had perforated 
fins and 3 did not. 

9. Core 5 produced results similar to 3 and 6 because of its 
extremely high friction factor even though the plate spacing is nearer 
to that of cores 1, 2 and 4. 

10. All cores produced a substantial increase in thermal effi- 
ciency over that obtainable with the open cycle. 
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Cj^=.1946 in 
C2=.1667 in 
Cj=.0431 in 
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FIGURE 11. 
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SECTIONAL DRAWING ILLUSTRATING CORE GEOMETRIES 
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Table I Characteristics and Program Results - Core 1 



Designation, Plain plate-fin surface 10.27 T 
Plate spacing, b « 0.544 in. 

Flow passage hydraulic diameter, 4rj, = .01259 ft. 
Compactness, jS - 289.93 ft^/ft^ 

Fin metal thickness = .010 in. 

Fin area/total area * 0.863 

Fin spacing at plates, c ® 0.1946 In. 



Results at thermal efficiency “ .34 



Frontal 


Length 


Regenerator 


Weight 




Area 


(ft) 


Effectiveness 


(lb) 


\ P / total 


3.0 

3.5 


desired 

6.2 


thermal efficiency 
.793 


not obtainable 
1122 


.116 


4.0 


4.4 


.754 


910 


.079 


4.5 


3.7 


.738 


816 


.062 


5.0 


3.3 


.731 


853 


.053 


5.5 


2.9 


.720 


825 


.046 


6.0 


2.7 


.719 


838 


.042 


6.5 


2.5 


.716 


840 


.038 


7.0 


2.3 


.710 


832 


.036 





Results 


at maximum thermal efficiency 






Frontal 


Thermal 


Optlmxim Regenerator 


Weight 




Area 


Efficiency 


Length (ft) Effectiveness 


(lb) 


\ P /total 



3.0 


.320 


5.9 


.765 


985 


.143 


3.5 


.343 


7.3 


.816 


1321 


.139 


4.0 


.359 


8.5 


.850 


1758 


.132 


4.5 


.371 


9.7 


.875 


2257 


.129 


5.0 


.382 


9.7 


.885 


2508 


.114 


5.5 


.391 


10.2 


.897 


2901 


.106 


6.0 


.399 


11.0 


.910 


3413 


.103 


6.5 


.406 


11.7 


.919 


3932 


.100 


7.0 


.411 


12.3 


.927 


4452 


.097 
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Table II Characteristics and Program Results - Core 2 



Designation, Plain plate-fin surface 12.00 T 
Plate spacing, b = 0.250 in. 

Flow passage hydraulic diameter, 4r^ = .009412 ft. 
Compactness, /S = 392.7 ft^/ft 
Fin metal thickness = .006 in. 

Fin area/total area =* .773 

Fin spacing at plates, c = .1667 in. 



Results at thermal efficiency = .34 



Frontal 


Length 


Regenerator 


Weight 






Area 


(ft) 


Effectiveness 


(lb) 


^ p ) 


total 


3.0 

3.5 


desired 

3.1 


thermal efficiency not 
.769 


obtainable 

607 


.095 




4.0 


2.6 


.749 


582 


.071 




4.5 


2.2 


.730 


554 


.057 




5.0 


2.0 


.7 26 


560 


.049 




5.5 


1.8 


.718 


554 


.044 




6.0 


1.7 


.719 


571 


.040 




6.5 


1.6 


.719 


582 


.037 




7.0 


1.5 


.718 


588 


.035 







Results 


at maximum 


thermal efficiency 






Frontal 


Thermal 


Optimum 


Regenerator 


Weight 


\ ^ /total 


Area 


Efficiency 


Length (f 


t) Effectiveness 


(lb) 


3.0 


.332 


4.0 


.803 


672 


. 149 


3.5 


.349 


4.9 


.838 


960 


.142 


4.0 


.365 


5.8 


.866 


1299 


.135 


4.5 


.377 


6.6 


.887 


1663 


.128 


5.0 


.388 


7.2 


.902 


2016 


. 121 


5.5 


.397 


7.8 


.914 


2402 


.116 


6.0 


.404 


8.3 


.924 


2788 


.111 


6.5 


.410 


8.3 


.928 


2984 


. 101 


7.0 


.415 


8.4 


.933 


3292 


.096 
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Table III Characteristics and Program Results -- Core 3 



Designation, Plain place-fin surface 46.45 T 
Plate spacing, b - 0.100 in. 

Flow passage hydraulic diameter, 4rj^ ~ 0.002643 ft 
Compactness, Q - 1332o45 fc^/fc^ 

Fin metal thickness ~ .002 in. 

Fin area/total area ~ .837 
Fin spacing at plates, c « .0431 



Results at thermal efficiency 



,34 



Frontal 

Area 



Length 

(ft) 



Regenerator 

Effectiveness 



Weight 

(lb) 



aP 
P 



total 



3.0 


desired 


thermal efficiency not 


obtainable 




3.5 


desired 


thermal efficiency not 


obtainable 




4.0 


0.5 


.819 


206 


.139 


4.5 


0.3 


.760 


139 


.084 


5.0 


0.3 


.779 


154 


.077 


5.5 


0.2 


,733 


113 


.055 


6.0 


0.2 


.750 


123 


.051 


6.5 


0.2 


.765 


134 


.048 


7,0 


0,2 


,778 


144 


.045 



Results at maximum thermal efficiency 



Frontal 


Thermal 


Optimum 


Regenerator 


Weight 


/A P ) 


Area 


Efficiency 


Length (ft) 


Effectiveness 


(lb) 




3.0 


.303 


0.5 


.769 


154 


. 186 


3.5 


.323 


0.6 


.820 


216 


.185 


4.0 


,342 


0.6 


.841 


247 


. 162 


4.5 


.355 


0.7 


.872 


324 


. 164 


5,0 


.368 


0.7 


.884 


360 


.148 


5.5 


.377 


0.8 


.905 


452 


.150 


6,0 


.389 


0.7 


.904 


432 


. 121 


6,5 


.397 


0.8 


.920 


535 


.124 


7.0 


.403 


0.9 


.933 


648 


.126 



total 
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Table IV Characteristics and Program Results - Core 4 



Designation, Plain plate-fin surface 6.2 
Plate spacing, b = .405 in. 

Flow passage hydraulic diameter, 4rj^ = 0.0182 ft 
Compactness, “ 204.0 ft2/ft^ 

Fin metal thickness = 0.010 in. 

Fin area/total area = .728 

Fin spacing at plates, c = .322 in. 



Results at thermal efficiency = .34 



Frontal 


Length 


Regenerator 


Weight 


/aP 


Area 


(ft) 


Effectiveness 


(lb) 


P 



3.0 


6.2 


.756 


887 


.082 


3.5 


5.8 


.736 


968 


.062 


4.0 


5.7 


.726 


1087 


.052 


4,5 


5.6 


.721 


1202 


.046 


5.0 


5.4 


.716 


1287 


.042 


5.5 


5.2 


.714 


1364 


,038 


6.0 


4.9 


.712 


1402 


.035 


6.5 


4.6 


.710 


1426 


.033 


7.0 


4.3 


.708 


1435 


,031 





Results 


at maximum thermal efficiency 




Frontal 


Thermal 


Optimum 


Regenerator 


Weight 


Area 


Efficiency 


Length (ft) 


Effectiveness 


(lb) 



3.0 


.355 


10.9 


.848 


1559 


3.5 


.371 


13,3 


.869 


2219 


4.0 


.382 


16.2 


.886 


3090 


4.5 


.391 


18.5 


,897 


3969 


5.0 


.399 


20.7 


.907 


4935 


5.5 


.405 


21.4 


.912 


5612 


6.0 


.411 


22.4 


.920 


6551 


6.5 


.415 


22.7 


.923 


7035 


7.0 


.420 


23.5 


.929 


7843 



total 



^ P 



total 



.137 
.123 
. 116 
.106 
.099 
.091 
.086 
.079 
.075 
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Table V Characteristics and Program Results - Core 5 



Designation, Wavy-fin plate~fin surface 11.5-3/8W 
Plate spacing, b = .375 in. 

Flow passage hydraulic diameter, 4r^^ - 0.00993 ft. 
Compactness , 6' = 345 ft 2/ft 3 
Fin metal thickness - .010 in. 

Fin area/total area - .822 

Fin spacing at plates, c - .087 in. 



Results at thermal efficiency ~ ,34 



Frontal 

Area 



Length Regenerator 

(ft) Effectiveness 



Weight 

(lb) 




total 



3.0 


desired 


3,5 


desired 


4.0 


desired 


4.5 


0.9 


5.0 


0.8 


5.5 


0.7 


6.0 


0.7 


6.5 


0.6 


7.0 


0.6 



themnal efficiency 
thermal efficiency 
thermal efficiency 
.111 
.758 
L739 
.745 
.721 
.726 



not obtainable 
not obtainable 
not obtainable 



335 


.092 


331 


.073 


319 


.060 


347 


.054 


323 


.046 


347 


.043 



Results at maximum thermal efficiency 



Frontal 

Area 



Thermal 

Efficiency 



Optimum Regenerator 

Length (ft) Effectiveness 



Weight 

Ob) 




total 



3.0 


.301 


0.9 


.742 


223 


. 168 


3.5 


.321 


1.1 


.788 


319 


.161 


4.0 


.338 


1.3 


.821 


430 


.152 


4.5 


.352 


1.6 


.854 


596 


,153 


5.0 


.364 


1.7 


,866 


703 


, 137 


5.5 


.374 


2.0 


.888 


910 


.131 


6.0 


.383 


2.2 


.900 


1092 


.131 


6.5 


.391 


2.3 


.907 


1237 


.120 


7.0 


.398 


2.5 


.916 


1448 


.116 



^Weights are low since they were based on core length not flow 
passage length 
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Table VI Characteristics and Program Results - Core 6 



Designation, Perforated-fin, plate-fin surface. Solar #2 
Plate spacing, b = 0.100 in. 

Flow passage hydraulic diameter, 4r, = .0046272 ft. 
Compactness, /O = 673.143 ft^/ft^ 

Fin metal thickness - .0047 in. 

Fin area/total area = .645 

Fin spacing at plates, c = .075 in. 



Results at thermal efficiency = .34 



Frontal Length Regenerator 

Area (ft) Effectiveness 



Weight /a P \ 

(lb) I p y 

total 



3,0 


desired 


thermal efficiency 


not 


obtainable 




3.5 


desired 


thermal efficiency 


not 


obtainable 




A.O 


desired 


thermal efficiency 


not 


obtainable 




4,5 


desired 


thermal efficiency 


not 


obtainable 




5.0 


0.7 


.768 




523 


.089 


5,5 


0 . 6 


,763 




493 


,074 


6.0 


0.5 


.749 




449 


.061 


6.5 


0.4 


.725 




389 


.051 


7.0 


0.4 


.741 




419 


.048 



Results at maximum thermal efficiency 



Frontal 


Thermal 


Optimum 


Regenerator 


Weight 


/aP'' 


Area 


Efficiency 


Length (ft) 


Effectiveness 


(lb) 


1 P 



3.0 


.287 


0,8 


.722 


359 


. 186 


3.5 


.305 


1.0 


,766 


523 


.179 


4.0 


.322 


1.1 


.793 


658 


. 164 


4.5 


.338 


1.2 


.822 


808 


.153 


5.0 


,352 


1.2 


,840 


897 


.136 


5.5 


.365 


1.4 


.871 


1151 


. 140 


6.0 


.375 


1.5 


.889 


1346 


.136 


6.5 


.386 


1.5 


.901 


1453 


, 125 


7.0 


.393 


1.6 


,913 


1675 


. 123 



total 
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LENGTH (FT) 




FRONTAL AREA (FT^) 

FIGURE 12. LENGTHS REQUIRED WITH VARIOUS FRONTAL AREAS TO 
OBTAIN 34 PERCENT THERMAL EFFICIENCY 
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CORE WEIGHT (LBS) 




2.0 3.0 4.0 5.0 6.0 7.0 

FRONTAL AREA (FT^) 



FIGURE 13. CORE WEIGHT VERSUS FRONTAL AREA FOR 34 PERCENT 
THERMAL EFFICIENCY 
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TOTAL RELATIVE PRESSURE DROP 




.69 .71 .73 .75 .77 .79 

REGENERATOR EFFECTIVENESS 

FIGURE 14. TOTAL RELATIVE PRESSURE DROP VERSUS REGENERATOR 
EFFECTIVENESS FOR -M .34 
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THERMAL EFFICIENCY 




2.0 3.0 4.0 5.0 6.0 7.0 

FRONTAL AREA (FT^) 

FIGURE 15. MAXIMUM THERMAL EFFICIENCY OBTAINABLE VERSUS 
FRONTAL AREA 
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LENGTH (FT) 




FRONTAL AREA (FT^) 



FIGURE 16 . LENGTHS REC3UIRED WITH VARIOUS FRONTAL AREAS FOR 
MAXIMUM THERMAL EFFICIENCY 
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CORE WEIGHT (LBS) 




2.0 3.0 4.0 5.0 6.0 7.0 

FRONTAL AIEA (FT^) 

FIGURE 17. CORE WEIGHT VERSUS FRONTAL AREA FOR MAXIMUM 
THERMAL EFFICIENCY 
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LEiMoTH = 1.0 REQUIRED FRONTAL AREA 
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SPECIFlL REYNOLDS REYNOLDS TRANSFER 

FUEL CON. NO. COLO NO. HOT UNITS CMIN/CMAX 



LENGTH = 4.3 REQUIKtD EkUNTAL AREA 




CL«| 


fsj 






0.-J 






< 


cn 




• 


< 


o 




<♦- 


o 


>0 




<K- 


o 




HO 


• 


H 




HO 


• 


H 


-JH 


o 




O 


-JH 


o 




LU 








UU 






Q 








o 







H 


• 












H 


• 




X 














X 


H 


its 


o 


O 


H 


>o 








0 


iTk 


H 


»— < 


o* 




^4 










a 




UJ 








X 






UJ 






3 








< 






JB 














X 


nO 


















o 


vO 
















• 




O' 










UJ 






o 


2: 


• 




UJ 






X 


nO 




o 


)-fl4 


0 




X 


LT^ 


NT 


o 


• 


H 


(N 


X 






o 


• 


H 


-J 








o 






J 


CO 




o 














0 


p<H 




> 














> 














a: 




















Lb 


















• 


LLCO 


«N 










{/) 






o 


00 »- 






to 






OC^ 




m 


r- 




• 




QLoO 




m 


OUJ 




H 




<2 


r\j 




OUJ 




H 


H-X 






^4 


QCZ> 






HZ 






<uu 








H 






<UJ 


>o 




QO 


m 












QC> 


0 




Lb*-' 














UJ^ 


H 






• 












ZH 


• 




'XIO 


o 






t/^ 






LbO 


0 




OUJ 






• 


QH 






OUJ 






LUUL 




(NJ 


ro 




• 


• 


LUUL 




CNi 


acILL 




H 


H 


OX 


m 


fO 


CLUL 




H 


UU 






X) 


X 


sO 




UJ 














> • 


CO 


II 
















LUO 


*-<< 


















OlX 




< 








a: 












UJ 


QC 


















a: 








<x 












<r 


< 






o 


O' 


m 


• 


oOQ 






0 


H 


m 




o 


CL 


ITS 








0 *-i 


0 


Qu 


H-h- 








JO 


• 


< 


HH 






< 


o 






OO 




H 


< 


0 




-Jo: 


• 








ro 


Z 


joe 


• 




Lb 


o 






> • 


0 


0 


Lb 


0 




o 








UJQ 


c\j 


LL 


0 






LL 






CVJ 


a:^ 




LL 


LL 
































• 






0 






> 4 * 






Q. 








UJ 






ca. 


> 






m 


• 




□c 


> 






o 








02: 




4 


O 






-JX 








— 0 


^;r 


o 


JZ 






<Uj 


o 






LLO 


0^ 


0 


<UJ 


0 




5 .*-* 


>r 






»— « 


ro 


Lb 




>T 




q:o 


m 






0 J 


• 


QC 


CLO 


fO 




UJ*— « 


• 




nO 


UJLU 


0 




Lb— • 


• 




XLL 


o 




• 


a-o 






XLL 


0 


m 


HLL 






00 


00 LL 






HLL 




CL 


•JJ 






lo 








LU 


















0 

• 




















lT> 






















JC 






o 






o 






M 


0 






-J 


• 




00 








J 


« 




LLUJ 




<N 


• 






X 


U.UJ 


0 


(NJ 


H 


'O 


a. 


00 






►— 


H 


H 


Cl 


uo <t 






m 






0 


LO< 


0 




oo: 












z 


00 CL 


(NJ 




< 


<M 










UJ 


< 


(NJ 
















-J 











Q. 










>v 










Q. J 


H 






• 


< 


0 




• 


CNJ 


<H 


^4 


0 


0 


0 


HO 


• 


H 


UN 


0 


JH 


0 




0 




UJ 










0 









• 














t 














H 


• 




0 






00 








X 


UN 


UN 


00 






>0 








0 


O' 


H 


>0 
















(J* 




^4 








X 






LU 








X 






< 






3: 








< 






X 


H 












X 


00 




0 


>o 












0 


>o 


• 


V. 












• 




( 7 ' 


0 


z 


• 




UJ 






0 


z 


• 


0 




0 




X 


(NJ 




0 


M 


0 


(NJ 


X 






0 


• 


H 


(NJ 


X 




CNI 


0 






J 


O' 




(NJ 


0 












0 


^4 


















> 














QC 














QC 






UJ 














UJ 




• 


ULlO 


0 










• 


LLlO 


vO 


•-J 


00 b- 






LO 






CO 


LO*— 


CVJ 




Z^ 


• 




OCoO 




ro 


0 


Zn-i 


• 




<z 


ro 




Olu 




H 


UN 


<z 


ro 




QCO 






l-Z 






^4 


QCO 






, H 






<UJ 


UN 






H 












QC> 


H 


















Lb*— < 


H 


















ZJ 


• 


















UJO 


0 






UO 




• 


OH 




UN 


(^UJ 






• 


OH 




rO 


JO 


• 


• 


LULL 




(NJ 


ro 


JO 


• 


H 


OX 


(NJ 


ro 


ocll 




H 


H 


OX 


vO 


00 


z 


0 




LU 






00 


z 


UN 




>- • 


c^ 


II 










> • 


^4 




LUO 














LUO 


(N 




QCZ 




< 










(XZ 










UJ 


QC 


















QC 








ers 






00 






< 


< 






0 






• 


LOO 






0 


UN 


UN 


• 


OOQ 




UN 


Q-J 




J 


o—t 


0 


CL 




o_i 




fH 


JO 


• 


< 


HH 


»-H 






JQ 


• 




00 


H 


»— 


<r 


0 






00 


H 




z 


0 


Z 


JQC 


• 






Z 


ro 




> • 




0 


Lb 


0 






> • 


ro 




LUO 


(N 


QC 


o 








LbO 


CNJ 


H 


QCZ 




LL 


LL 






(NJ 


QCZ 
























• 






CQ 








• 












UJ 






CL 


H 






UN 


• 




QC 


> 






UN 


• 






02 




*-ii 4 


0 








oz 






-.0 




0 


JZ 








*-0 


^0 




LLO 




Of 


<UJ 


00 






LLO 


CT' 




t — 1 


rO 


OJ 




ro 








ro 


UN 


0 J 


• 


X 


a:^ 


CO 




(O' 


0 J 


• 




LULU 


0 




UJ — 


• 




UN 


LUUJ 


0 


• 


CLO 






XU. 


0 


ro 


• 


a.3 




00 


LOX 






HLL 




CL 


00 


LOLL 




m 








Lb 






UN 












u> 

• 




















UN 






















JC 












O' 






M 


o 






0 






00 








J 


• 




00 






• 






X 


LLlU 


O' 


(NJ 


• 






00 






H 


H 


QO 


Cl 


QU 






UN 






0 


U0< 


UN 




UN 












Z 


00 QC 


(NJ 
















LU 


< 


















J 















75 



LENGTH = 6.0 RhgUIRED FRONTAL AREA 
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DESIRED THERMAL EFFICIENCY NOT OBTAINABLE WITH FRONTAL AREA 
MAXIMUM THERMAL EFFICIENCY =.320 
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393 1975. 1808. 2.90 0.965 



FRONTAL AREA = 4.5 REQUIRED LENGTH 
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APPENDIX A 



Computer Program Usage 

The computer program is written in FORTRAN IV language described in 
(9). The main program and all subroutines were written using double 
precision in order to obtain the best possible accuracy. This is parti- 
cularly essential in subroutines where least squares curve fitting is 
utilized to determine a functional relationship from a set of data points. 

Required computer storage capacity is small as no subscripted vari- 
ables or arrays were introduced while writing the program. Essentially, 
the only storage requirement is for program compilation. 

Approximately three minutes of computer time on an IBM System/360 
is required for each output type described in section six with the excep- 
tion of output types one and six which require only about 30 seconds. 

a. Subroutine XMAIN 

The main program is designed for sizing counterflow heat 
exchangers only. Subroutine XMAIN is general however and could be uti- 
lized for numerous other flow arrangements provided that function GENFX 
(TU, RATIO), called by subroutine XMAIN, is changed to a regenerator ef- 
fectiveness function corresponding to the desired flow arrangement. 

b. Subroutines Dependent on Internal Core Geometry 

The relationships of friction factor and Colburn j -factor 
vs. Reynolds number must correspond to the core geometry being investi- 
gated. These relationships are fed into the program as function sub- 
programs FRICC(RE), FRICH(RE), STPRC(RE) AND STPRH(RE) . An example is 
shown on page 111. Only the cold side functions are shown as the same 
core was used for both hot and cold sides in all cases. 

Subroutine KEKC (SIGMA, RE, XKE, XKC) , and as a function of 
(X“ and Re, depends on the shape of the core flow passages. Subroutines 
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were written for square and triangular passages and are shown on pages 
105 and 106. 

c. Data Cards 

All general inputs for output types one through five 
are read into the computer program from the first four data cards. Ten 
columns are allotted for each input. Each of the general inputs must 
be of floating point form. The inputs with corresponding units are 
listed below in the proper order: 

compressor inlet temperature ®R 

-2 

compressor inlet pressure Ibf-ft 

compressor pressure ratio 

compressor efficiency 

turbine inlet temperature 

turbine efficiency 

assumed regenerator effectiveness 

fuel inlet temperature ®R 

combustor efficiency 

lower heating value of fuel BTU-lbm 

hot side regenerator exhaust pressure lbf~ft 

horsepower output hp 

assumed relative pressure drop 
through regenerator cold side 

assumed relative pressure drop 
through regenerator hot side 

gas constant for air 

regenerator plate thickness 

cold side transfer area/volume 
between cold side plates 



-1 



ft-lbf(lbm-^R) 



-1 



in 



ft 



87 



hot side transfer area/volume 
between hot side plates 

cold side plate spacing 

hot side plate spacing 

cold side hydraulic radius 

hot side hydraulic radius 

mass flow rate 

relative pressure drop through 
combustor 

cold side fin thickness 

hot side fin thickness 

fin thermal conductance 

cold side fin or extended area/ 
total area 

hot side fin or extended area/ 
total area 



ft-lbf(lbm-°R)‘^ 



in 

in 

ft 

ft 



Ibm-hr 



-1 



ft 

ft 

BTU-(hr-ft^-'’F/ft) 



-1 



-3 



density of regenerator material Ibm-ft 

If output type six is desired, only the first 12 general inputs are 
necessary. The third and fourth data cards must be placed in the deck; 
however, they may be left completely blank. The data cards utilized in 
generating the results for core 2 are shown on page 95 . 

The fifth data card is used to designate the type of output desired. 
The fixed point integer "1** in column one designates output type one; 
the fixed point integer "2** in column one designates output type two 
etc. 



The sixth and last data card depends on the type of output designat- 
ed. In all cases numbers must be in floating point format. Ten columns 

are allowed for each input beginning in column one. Units are either ft. 
2 

or ft. corresponding to length or area. The inputs for each type out- 
put are listed below in the proper order: 
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output type two; 



output type three; 



output type four: 
output type five: 



output type one: length 

frontal area 
shortest length 
longest length 
length increment size 
thermal efficiency 
smallest area 
largest area 
area increment size 
thermal efficiency 
thermal efficiency 
smallest area 
largest area 
area increment size 
no data card 
d. Incidental Instruction in Program Usage 

There are several potential program trouble areas that 
the user should be aware of when investigating various core geometries. 
First of all, it is important to keep in mind the range of Reynolds 
numbers over which the functional relationships of friction factor and 
Colburn- j factor have been obtained. Frontal areas must be chosen so that 
the hot and cold Reynolds numbers fall in this range. This is of course 
largely dependent on internal regenerator core geometry and mass flow 
rate. Secondly, experience has shown that the iteration process diverges 
when the total relative pressure drop through both sides of the core and 
combustor exceeds approximately 307o* This pressure drop is mainly a 
function of regenerator length and core plate spacing. In general, shorter 



output type six: 
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lengths are necessary as plate spacing is decreased* 

Main program statements 061, 125, 126 and 141 designate the initial 
frontal area and length used in incrementing out the desired values for 
output types two and four. These statements, as shown in the program 
printout, are set for a 600 hp gas turbine with a regenerator having a 
core plate spacing of .2 inches or greater. For cores with plate spacing 
smaller than ,2 inches or exceptionally high friction factor it may be 
necessary to shorten the length in statement 125 and to increase the 
frontal area in statements 061, 126 and 141 to get the incrementing pro- 
cess started properly It may also be necessary to increase the frontal 
area in the above statements for larger gas turbine arrangements* 

Main program statements 62, 92, 125, 140, 153 and 177 designate the 

length and area increments to be used in the main program output types. 

2 

They are presently set for .1 ft and .1 ft but may be changed as neces- 
sary if better accuracy is desired. 



e. Listing of Program Variables 



A1 



dummy variable used in first read statement 



A2 



dummy variable used in first read statement 



AFINC cold side fin area/total area 



AFINH hot side fin area/total area 



AFLOWC cold side free-flow area 



AFLOWH hot side free-flow area 



AFRINC incremental frontal area 



AFRON 



frontal area 



AFRONl smallest frontal area considered 



AFR0N2 largest frontal area considered 



AHTC 



cold side heat transfer area 



AHTH 



hot side heat transfer area 



90 



ALFAC 


total transfer area of cold side of 
regenerator/ total regenerator volume 


ALFAH 


total transfer area of hot side of 
regenerator/total regenerator volume 


BC 


cold side plate spacing 


BETAC 


cold side transfer area/cold side volume 


BETAH 


hot side transfer area/hot side volume 


BH 


hot side plate spacing 


CAPO 


cold side capacity rate 


CAPH 


hot side capacity rate 


CEFF 


compressor efficiency 


COBJC 


cold side j factor 


COBJH 


hot side j factor 


COMEFF 


combustor efficiency 


CPC 


specific heat of cold side air 


CPH 


specific heat of hot side air 


CRATIO 


compression ratio 


DEN 


density of core material 


FC 


cold side friction factor 


FCOND 


fin thermal conductance 


FH 


hot side friction factor 


FINEFC 


cold side fin efficiency 


FINEFH 


hot side fin efficiency 


FINTHC 


cold side fin thickness 


FINTHH 


hot side fin thickness 


GC 


cold side mass flow rate per unit area 


GENEFF 


regenerator effectiveness 


GH 


hot side mass flow rate per unit area 


HI 


enthalpy of air at compressor inlet 
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H2 


compressor discharge enthalpy 


H3 


combustor inlet enthalpy 


H4 


turbine inlet enthalpy 


H5 


turbine exhaust enthalpy 


H2S 


compressor discharge ideal enthalpy 


H5S 


turbine exhaust ideal enthalpy 


HEATC 


cold side heat transfer coefficient 


HEATH 


hot side heat transfer coefficient 


HETVAL 


lower heating value of fuel 


HF 


fuel enthalpy 


HP 


horsepower output 


PI 


compressor inlet pressure 


P2 


compressor discharge pressure 


P3 


combustor inlet pressure 


P4 


turbine inlet pressure 


P5 


turbine exhaust pressure 


P6 


regenerator hot side exhaust pressure 


P23 


pressure drop through regenerator cold side 


P56 


pressure drop through regenerator hot side 


P23R 


relative pressure drop through regenerator 
cold side 


P34R 


relative pressure drop through combustor 


P56R 


relative pressure drop through regenerator 
hot side 


PRl 


relative pressure at compressor inlet 


PR2S 


relative pressure at point 2s. 


PR5S 


relative pressure at point 5s. 


PR4 


relative pressure at turbine inlet 


PRC 


cold side Prandtl number 
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PRH 


hot side Prandtl number 


R 


gas constant for air 


RADC 


cold side hydraulic radius 


RADH 


hot side hydraulic radius 


RATIO 


capacity rate ratio 


REG 


cold side Reynolds number 


REH 


hot side Reynolds number 


RELP 


total relative pressure drop 


RFA 


fuel to air ratio 


SFC 


specific fuel consumption 


SIGMAC 


cold side free flow area/cold side 
frontal area 


SIG3MAH 


hot side free flow area/hot side 
frontal area 


STC 


cold side Stanton number 


STH 


hot side Stanton number 


SUREFC 


cold side surface efficiency 


SUREFH 


hot side surface efficiency 


SVOL 


solid core volume 


T1 


compressor inlet temperature 


T2 


compressor discharge temperature 


T3 


combustor inlet temperature 


T4 


turbine inlet temperature 


T5 


turbine exhaust temperature 


T6 


regenerator hot side exhaust temperature 


T2S 


temperature at point 


T5S 


temperature at point 


TCM 


mean cold side temperature 


TEFF 


turbine efficiency 
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TF 

THEFF 

THM 

TU 

U 

V2 

V3 

V5 

V6 

VCM 

VHM 

VIC 

VIH 

VOL 

W 

WT 

XKCC 

XKCH 

XKEC 

XKEH 

XL 

XLl 

XL2 

XLINC 

XMDOT 



fuel temperature 

thermal efficiency 

mean hot side temperature 

number of heat transfer units 

overall heat transfer coefficient 

cold side inlet specific volume 

cold side outlet specific volume 

hot side inlet specific volume 

hot side outlet specific volume 

mean cold side specific volume 

mean hot side specific volume 

cold side viscosity 

hot side viscosity 

regenerator volume 

regenerator plate thickness 

regenerator weight 

cold side contraction coefficient 

hot side contraction coefficient 

cold side expansion coefficient 

hot side expansion coefficient 

% 

regenerator length 
shortest length considered 
longest length considered 
length increment 
mass flow rate 
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APPENDIX B 



Curve Fitting 

It is frequently necessary to represent by a functional relation- 
ship data that is given in the form of a set of x-y points. One may 
want to know the value of "y" corresponding to some "x" that lies between 
two tabulated "x" values or may want to find the "y” value corresponding 
to an "x" value outside the range of the tabulated "x" values. The 
problem then is to find a curve that approximates the given data with 
sufficient accuracy. 

Before the computer program could be fully developed, various func- 
tional relationships of the type mentioned above had to be obtained. 

Many of the necessary relationships were available in literature? how- 
ever, seven of them were not and had to be obtained by some appropriate 
curve fitting method. The "least squares" method was utilized in fitting 
six of the curves and straight line approximations were used on the re- 
maining one. It was decided to use the least squares method because of 
its simplicity and the smoothness of fit that is realizable in most 
cases. 

a. Least Squares Theory. 

The deviation at a data point is defined as the difference 

between the tabulated y-value and the y-value computed from a functional 

relationship. The least squares method asks for a minimum of 

N 

(B-1) 

(°=l 

where N equals the number of data points, y' denotes a value computed 
from the functional relationship and y denotes the actual tabulated value. 
This leads to equations that in many oeses of praotieal intorest are 
linear and, in pfinctple at least easy to solve. If It ii assumed that 



the functional relationship desired is of the form 

a» + a,x + ^ (b-2) 



a , a-, a, must be determined to minimize 

0 1/ k ^ 

S=Y. (b-3) 

L.'^\ 

Taking the derivative of S with respect to and setting it equal to zero 



Lc^\ ^ 

a/ 2, 

- ^ {<^0 + >^L -\r O^lXc ‘ 0.K Xe — 

C-l 

N 

~ 2. H-£^, Xt’ "V X-i* -t~ ' * ' 

a . ' 



/\/ 



6-/ 



N •■ 2. AT -— ’ 

— ^ <51,9 ^z.Xt” -f- ' ' ' ■\^^<X(f — —O 

c-\ c-i r^i L^i 



A/ 

I' 

I'-i 



A/ ^ A/ 



or 



X/Ao +■ 



^ ( I ^ 1.1 I ^ \ ^ 

+ y^ Xt“ I ^1. 1 z I 

L‘-l ' ' ' /T; 



(B-4) 



This gives one equation. By taking the derivatives of S with respect to 



the other coefficients, a^ through a^^, another "k" equations are 





N 

i^t 



(B-5) 





fJ 




0 

L 
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The k+1 equations can then be solved for the k+1 coefficients. 

The value selected for k can be any value less than N* For best 
results, however, k should be made only large enough to give desired 
accuracy. As k is increased the computer round-off errors inherent in 
solving the set of simultaneous equations becomes more significant and 
accuracy is lost not gained. The generated curve may also tend to fluc- 
tuate between data points if k is too large, making a smooth fit impos- 
sible . 

b. Gauss Elimination. 

The method of Gauss elimination was used for solving the 
simultaneous equations mentioned above. A system of three equations in 
three unknowns will be used to illustrate the procedure because the 
nomenclature becomes somewhat complicated for the general case of **n** 
equations and **n** unknowns. Let and be the unknowns and the 

equations are 

(X\\X\ 13 ~ 

(B-6) 

The equations must be arranged so that a^^ is not zero. Next define a 
multiplier. 



/Mz- 



1 

^ II 



(B-7) 



The first equation is then multiplied by m 2 , and subtracted from the 
second equation. The result is 
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(B-S) 


But 

(57T ^ ^ 

so Xj^ has been eliminated from the second equation. 
Define : 


(B-9) 


/ 


(B-10) 


A. “ A ^ 


(B-11) 


then equation (B-8) becomes 

"d" ^>3 ^ A. 

If a multiplier is now defined for the third equation. 


(B-12) 


^3 1 


(B-13) 



and the first equation is subtracted frotii the third, the result is 



^3Z -\*^33^Xj ■= 

where 


(B-14) 




(B-15) 


(A 32 ' ~ <^S3 ~ ‘^3 O- 13 


(B-16) 




(B-17) 
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The resulting three equations in three unknowns are 



X| -h 13 Xj — 

^ i 2.^ ^ 2. 'irj, 

(^%X Xz -{- Xj •= £.J 



If is now eliminated from the last equation using the same procedure 
the following is obtained: 

(k\{ X| + + 

ai-L'Xzf (B-19) 

X-s = 

It is now a straight forward process to solve for x^, to substitute that 
result into the second equation to find and finally to substitute in- 
to the first equation to find Xj^. 



0 >! 



"v - f ^v! 





•'0-[xX\ ^^13 X 3 ) 
(X.W 



(B-20) 



c. Curve Fitting Computer Program 

The computer program shown on the following pages was 
utilized in fitting the necessary curves. The main program manipulates 
the data point information into the equation form shown in part "a" of 
this appendix. Subroutine GAUSS 2 is then called and the equations are 
solved simultaneously using the method described above (10), The coeffi- 
cients, determined for a polynomial of the form y “ a^^ + 



®2^ + 33^ 



k-1 



aj^x are printed out along with a column of x-y 
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points generated from the polynomial. These points can then be compared 
with the original data points to detemine whether or not the fit is 
sufficiently accurate. 

A maximum of fifty data points can be read into the computer. For 
best results, they should be as evenly spaced along the "x" axis as pos- 
sible. A special problem arises when using the program to fit Colburn- j 
and friction factor vs. Reynolds number data. This data is normally 
plotted on logarithmic coordinates and is therefore not evenly spaced on 
linear coordinates. To remedy this situation, the natural logarithms of 
the "x" and "y" values can be taken before starting the curve fitting 

procedure. The resulting polynomial is of the form ln(y) - a^^ + a 2 (ln x) + 
2 k-1 

a 2 (ln x) logarithms can then be taken to deter- 

mine the desired "y” value. These additional steps are incorporated into 
the computer program as it is shown. If the user desires to fit data 
plotted on linear coordinates, cards 150, 160, 510, and 520 must be re- 
moved and card 530 must be changed to read: 

WRITE (6,55) XI (LQ), Yl(LQ) 

Coefficients shown in parts 1 and 2 of Table VII are for polynomials 
of the first form shown above. The remaining coefficients in Table VII 
are for polynomials of the form just described. 

When utilizing the program, the following data cards are necessary: 

(1) The number of data points to be used in fitting the 
curve is punched into the first three columns on the first data card. 

The order of the polynomial desired is punched into the second three 
columns of the first data card. Both of these numbers must be in fixed 
point integer form. The order of the polynomial can be any number up to 
one less than the number of data points. 
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(2) The x-y values of the data points are punched into 
the next N data cards, one data point per card. The x-value is punched 
into the first ten columns; the y-value is punched into the second ten 
columns. Both numbers must be of floating point form. 

(3) As mentioned above, the printout includes a column 
of x-y points generated from the computed polynomial. The first point 
is computed using the smallest x-value and the last point is computed 
using the largest x-value of the original data points. The size of the 
desired increments between these two data points is punched into the next 
data card in the first ten columns. Floating point form must be used. 

(4) If the user desires to fit more than one set of 
points he can place the next set of data cards, for those points, be- 
hind the set just described. 

(5) The next and last data card must contain the integer 
2 punched into column three. This signals that all data point sets have 
been fit and stops the program. 

d. Curve Fitting Results 

A sample of least squares curve fitting is shown in Figure 

18. All data points used in curve fitting and the resulting polynomial 
coefficients are shown in Table VII. 

The curve fit by straight line approximations is shown in Figure 

19, with the straight lines drawn in. Two straight lines were used, one 

covering the range,. 3 0*" .55 the other covering the range 

.55<^0~<C.75. This approximation was used rather than the least squares 
theory because small errors in cause virtually no error in the comput- 
er program results. The curve is from Kays (5). 



/ < 
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Table VII Data Points Used in Curve Fitting and Resulting 

Polynoi'nial Coerficlents Used In Prograi.'! Subroutines 



1 . 



T^R 



DynaiTiic Viscosity of Air ac Atmospheric Pressure from 
Hilsenrath (6 k )/^ ~ 17 i6 x i poise 

y/ 

Polynomial Coefficients 



180 


0.403S 


A] 




- 132588012 


D-01’ 


540 


1.7742 






. 2706893585 


n-02 


900 


1.556 


a3 




-.1318733093 


D-05 


1260 


1.7 58 






.7538864121 


D-12 


1620 


2,271 


A 5 




"^21455415 


D-12 


1980 


2 . 562 


a5 




^7 52 18 7450 


D-15 


2340 


2.824 


a7 




. 2117993215 


D-18 


2 700 


3.066 


a8 


- 


-.3862860359 


D-22 


3060 


3.290 


A9 




, 2782105119 


D-26 


■■'>20 


3,501 


A IC 




. 67b3288488 


0-32 



2, 


Prandtl Number of Air f 


rom Kiisenrath (6 ) v 




T°R 


Prandtl No. 


Polynomial Coefficients 


180 


0.770 


Al = 


. 7865630968 


0-00 


350 


0.739 


A2 = 


-.28/1584333 


D-04 


540 


0. 708 


A3 = 


- . 3666074673 


D-06 


720 


0.689 


A4 == 


-.1373281188 


D-09 


900 


0.680 


A5 = 


.1569904076 


D-Il 


1080 


0.680 


A6 - 


-. 1973361281 


D-14 


1260 


0.684 


A7 = 


.1041204710 


D-17 


1440 


0.689 


A8 = 


-2047252924 


D-21 


1620 


0.696 


A9 = 


-.1538501631 


D-25 


1800 


0.702 


AlO = 


.7960530427 


D-29 


3. 


Core- 1 (Plain plate-fi'> 


surface 


10.27 I) Friction 



from Xays and London {1 



Reynolds 

Number 


Friction 

Factor_ 


ll. 


lynomial Coef tie lent 


300 


.066 


Ai 


= .1060^^2201 


D+03 


400 


.051 


Al 


= -.7242808930 


D-K)2 


500 


.042 


A3 


- .1760821881 


D-H32 


700 


.032 




= ••.1171044519 


D+Ol 


900 


.026 


A5 


= -.2766351384 


D+00 


1250 


.0185 


A6 


-- .6028901137 


D-01 


1500 


0165 


A7 


- -.1810176966 


D-02 


2000 


.0135 


a8 


- - 7728780735 


D-03 


3000 


.011 


A9 


^ ,1143227913 


D-03 


4000 


.0095 


AiO 


= -.6409201873 


D-05 


6000 


.0082 


All 


= .1272474507 


D-06 



A1 = -.4132588012 x 10'^ 
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Core 1 (Plain plate-fin surface 10.27 T) Coburn- j 
Factor from Kays and London (1) 



Reynolds 

Number 



Coburn- j 
Factor 



Polynomial Coefficie a t s 



300 


.016 


A1 = 


-.9478036457 


D+01 


400 


,013 


A2 = 


.5133725600 


D+01 


500 


.0115 


A3 = 


-.1236662988 


D+Ol 


700 


.0087 


A4 = 


-.2136908795 


D+00 


900 


,0072 


A5 = 


. 1534935252 


D+00 


1250 


.006 


A6 = 


-. 2255QQ'='^<^0 


n-oi 


1500 


.00515 


A7 = 


-.3614040571 


D-03 


2000 


.0045 


A8 = 


.3150813943 


D-03 


3000 


.00375 


A9 = 


-.1222790555 


D-04 


4000 


,00345 


AlO = 


-.1473493662 


D-05 


6000 


.0032 


All = 


.9202049786 


D-07 


5, 


Core 2 (Plain plate-fin surface 12,00 T) 


Friction 




Factor from Kays and London 


(1> 




Reynolds 


Friction 








Number 


Factor 


Polynomial Coefficients 


200 


,080 


A1 = 


.2855654237 


D'K)3 


300 


.057 


A2 = 


-.4366648607 


D+02 


500 


.037 


A3 = 


-.7695241634 


D+02 


700 


.027 


A4 = 


.3558930941 


D-H)2 


1000 


.019 


A5 = 


-.5846569710 


D-fOl 


1500 


.015 


A6 = 


.2516782204 


D+00 


2000 


.013 


A7 = 


.3534056422 


D-01 


3000 


.0115 


A8 = 


-.3719622996 


D-02 


4000 


.0112 


A9 = 


-.3919235358 


D-04 


6000 


.0095 


AlO = 


.1641370472 


D-04 


8000 


.0036 


All = 


-.5085864719 


D-06 


6. 


Core 2 (Plain 


plate-fin surface 12.00 T) 


Coburn- j 




from Kays and 


London ( 1) 






Reynolds 


Coburn- j 








Number 


Factor 


Polynomial Coefficients 


200 


,016 


A1 = 


-.1678025363 


D+04 


300 


,012 


A2 = 


. 1147089449 


D-fOA 


500 


,008 


A3 = 


-.2302525106 


D+03 


700 


,0064 


A4 = 


-.9930779398 


D+01 


1000 


,0049 


A5 = 


.8790472874 


D-l-01 


1500 


,0038 


A6 = 


-.7254951992 


D+00 


2000 


.00355 


A7 = 


-.5905850854 


D-01 


3000 


.0036 


A8 = 


.6125863884 


D-02 


4000 


,0035 


A9 = 


.9700663367 


D-03 


6000 


,00325 


AlO = 


-.1394667289 


D-03 


8000 


.003 


All = 


.4807945572 


D-05 



Factor 
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Core 3 (Plain piate-fin surface 46.45 T) Friction Factor 
from Kays and London i - 



Pveyaolds 

Number 


Friction 

Factor 


Folynomi al Coe f f ici ent s 


300 


.059 


A1 


St 


.5218534745 


I>-K)2 


400 


.045 


A2 




-.3337432034 


D+02 


500 


.037 


-A3 




,8259113372 


D+Ol 


700 


.0265 


A4 




-.1802267981 


D+Ol 


900 


.020 


A5 




.4388704605 


D-K)0 


1250 


.0163 


a6 


= 


-.4819876576 


D-Oi 


1500 


.0135 


A7 


- 


-.3056062581 


D-02 


2000 


.0128 


AS 




. 7437300503 


D-03 


3000 


.0115 


A9 




.3659712646 


D-04 


4000 


.0112 


A 10 




-.1134921023 


D-04 


uOOO 


.0110 


All 


St 


.5000930552 


D-0.5 



3. Core 3 {Plain plate-fin surface 46,45 T) Coburn- j Factor 
from Kays and London (i;“ 



Reynolds Cobum- j 



>3 umber 


Factor 


Po 


iynomial Coefficients 


300 


.0135 


A1 


a 


-.9179528233 


!>l-01 


400 


,0105 


A 2 




.2549082428 


D+01 


500 


.0036 


A3 




-.1138319133 


D-fOl 


700 


.0065 


A4 


- 


. 1312317457 


D+01 


900 


.0053 


A5 




-.5200720093 


D+00 


1250 


.0041 


A6 


Si 


.6707637499 


D-01 


1500 


.0036 


A 7 




.2230241152 


D-02 


2000 


.00295 


A8 


JS 


-.1010311829 


D-02 


3000 


.0026 


A9 




. 1206538219 


D-04 


4000 


.0025 


A 10 


=: 


.8449625214 


D-05 


6000 


.00225 


All 


Si 


-.4529482542 


D-06 



9. Core 4 (Plain plate-fi' surface 6.2) Friction Factor 
from Kays and London «1 



Reynolds 

Number 


Friction 

Factor 


600 


.027 


800 


.022 


1250 


.015 


2000 


.010 


2300 


.0098 


2700 


.0093 


3000 


. 0092 


4000 


.00875 


6000 


,0081 


8000 


.00775 


10000 


.0073 



Folynomlal Coef ficients 



Ai 


IT 


,4512980279 


D+05 


A 2 


- 


-.3083071127 


D+05 


A3 


Si 


.7727049904 


D+04 


A4 


S. 


-.7410141046 


D+03 


A5 




-.6317170678 


D+01 


A6 




,388064?509 


D+01 


A? 




.6306627426 


D-01 


A8 


ss 


.3219487556 


D-01 


A9 




-.1290508496 


D-01 


AlO 




. 1130983538 


D-02 


All 


s= 


-.3189159559 


D-04 



1.23 



10. Core 4 (Plain plate-fin surface 6.2) Coburn- j Factor 
from Kays and London (1 



Reynolds 

Number 



Coburn- j 
Factor 



Polynomial Coefficients 



600 


.0072 




A1 


= .7203819966 


D-H35 


800 


.0059 




A 2 


= -.4396069217 


D-H35 


1250 


.0043 




A3 


= .1220436223 


D+05 


2000 


,0032 




A4 


= -.1162973846 


D+04 


2300 


.003 




A5 


= -.1005649404 


Dt-02 


2700 


,0031 




A6 


= ,6019720585 


D+01 


3000 


.0033 




A7 


= .9239757582 


D-01 


4000 


.0034 




A8 


= ,5306752808 


D-01 


6000 


.0033 




A9 


= -.2033142744 


D-01 


8000 


.00325 




AlO 


= .1762226828 


D-02 


10000 


.00315 




All 


= -.4934822700 


D-04 


11. 


Core 5 (Wavy-fin plate-fin su-^face 11.5-3/8W) Fr: 




Factor from Kays 


and 


London 


(1 




Reynolds 


Friction 










Number 


Factor 




Polynomial Coefficients 


400 


,14 




A1 


= -.1703315790 


D^-05 


500 


.125 




A2 


= .1290411355 


D+05 


700 


.098 




A3 


= -.3694551557 


D-K)4 


900 


.083 




A4 


= .4359291610 


D+03 


1250 


.078 




A5 


= -.1981584997 


D+01 


1500 


.070 




A6 


= -.3594302044 


D+01 


2000 


,0625 




A7 


= .7340606854 


D-01 


3000 


.0525 




A8 


= .3091640055 


D-01 


4000 


.0475 




A9 


= -.1603543052 


D-02 


5000 


.043 




AlO 


= -.7179253410 


D-04 


7000 


,0375 




All 


= .5035665432 


D-05 


12. 


Core 5 (Wavy-fin 


plate-fin su'^face 11.5-3/81/) Col 




Factor from Kays 


and 


London 


(1) 




R.eynolds 


Coburn- j 










Number 


Factor 




Polynomial Coefficients 


400 


.019 




A1 


= ,1845757550 


D+05 


500 


,018 




A2 


= -.1406316423 


D-^05 


700 


.0175 




A3 


= 4071246982 


D+04 


900 


.01625 




A4 


= -.4950898847 


D+03 


1250 


,015 




A5 


= .5485777432 


D+01 


1500 


,014 




A6 


= .3828321238 


D+01 


2000 


.013 




A7 


= -.1208492523 


D+00 


3000 


.012 




A8 


= -.2934686653 


D-01 


4000 


.0096 




A9 


= ,1796872494 


D-02 


5000 


.009 




AlO 


= .5114973685 


D-04 


7000 


.0079 




All 


= -.4408556059 


D-05 
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13, Core 6 (Perforated -fin plate fin surface, Solar .- 4 ,) 
Friction Factor fro.fi Prcci . cl.fS). 



Reynolds 
\ e t 


Friction 

Factor 


150 


, 15 


200 


. 11 


300 


,073 


5G0 


.049 


700 


.035 


900 


.029 


1200 


.0225 


1500 


.020 


2000 


.020 


2500 


.020 


3C0O 


.020 



14, 


Core 6 (Perf orated -f 
Coburn-j Factor from 


Reynolds 


Coburn- j 


Number 


Factor 


150 


.04 


200 


.03 


300 


.019 


500 


.0115 


700 


.008 


90C 


.0065 


1200 


.0054 


1500 


. 0054 


2000 


.0065 


2500 


. 0063 


3000 


. 0058 



Polynomial Coe f f ic i eats 



Al 




.4120158900 


DH)4 


A 2 




-.3513598265 


D-K)4 






, 1039 787 565 


J-^04 


A4 


= 


-.32c32840240 


D-!'02 


A5 




-.2504768472 


D402 


a6 




.6576122775 


D+0 1 


A? 




-.5125072877 


D-lOO 






-.7245182215 


D-02 


A9 


* 


.3054392764 


D-02 


A 10 


sr 


1198421133 


D-03 


All 


= 


-.3775033884 


D-0? 



n plate fin surface, Solar #2) 
Pucci et . ai , ( 8 ) , 



Polynomial Coefficients 



Al 




-.7987254165 


D+04 


A2 




.6330003623 


D-K34 


A3 




-.1701306104 


D+04 


A4 


=. 


.9227812576 


D+02 


A5 




.3735667832 


D+02 


a6 


- 


-.6449659103 


D+01 


A7 




.1419750320 


D+00 


a8 




.2445179921 


D-01 


A9 


- 


.2010489258 


D-02 


AlO 


:=r 


-.5001272149 


D-03 


All 




.2068223608 


D-04 



PRANDTL NUMBER OF AIR 




FIGURE 18. DATA POINTS FIT .WITH POLYNOMIAL USING 
LEAST SQUARES THEORY 
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FIGURE 19. C VERSUS <T“FIT WITH STRAIGHT LINE APPROXIMATIONS 
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